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PARAMETRIC STUDY OF A FRANGIBLE-TUBE ENERGY-ABSORPTION 
SYSTEM FOR PROTECTION OF A NUCLEAR AIRCRAFT REACTOR 
by Richard L. Puthoff and Klaus H. Gumto 
Lewis Research Center 

SUMMARY 

The reactor containment vessel of a nuclear airplane must contain the fission prod- 
ucts even at ground impact. In this study frangible-tube energy absorbers are used to 
absorb the energy at impact and to prevent rupture of the containment vessel. Frangible 
tubes were selected because tests show that they absorb more energy per pound of ab- 
sorber than other known energy absorbers. 

A parametric study for determining the minimum system weight was made of an 
omnidirectional energy-absorption system utilizing frangible tubes. The system, com- 
posed of tubes, dies, and supporting structure, typified that which might be applied to a 
nuclear airplane. The results of the parametric study are the following: 

1. Increasing the absorber angular capability or increasing the specific energy of 
a frangible tube substantially decreases the energy absorber system weight. Variables 
of lesser effect but still of significance are the frontal impact velocity and the decelera- 
tion of the protected package. 

2. Changing the number of tubes or the ratio of tube diameter to die forming radius 
has a minor effect on system weight. 

From the parametric curves generated, a frangible-tube energy-absorption system 
would weigh 50 percent of the powerplant weight. This assumes an impact velocity of 
400 feet per second (122 m/sec) in the frontal direction and 250 feet per second 
(76 m/sec) in all other directions. 

The most promising method to reduce weight appears to be increasing the energy 
absorber operating angle and design integration of the reactor, shielding, containment 
vessel energy absorber, and aircraft structure. In this manner the system component 
serves both in its intended use during normal system operation and as an energy ab- 
sorber during impact. 



INTRODUCTION 


Nuclear airplanes have the potential to provide longer range than chemical aircraft. 
However, before the nuclear airplane can be used it must be accepted as being safe. 

Some of the safety problems connected with a nuclear aircraft are similar to a 
ground-based powerplant while others are uniquely different. The major difference is 
that the aircraft reactor is in motion and has the potential of crashing. To permit oper- 
ation near populated areas the containment vessel must remain intact and contain all fis- 
sion products upon impact with the ground or another aircraft. The containment vessel 
and the reactor contained within it have an enormous amount of kinetic energy when 
traveling at high speeds. To prevent rupture of the containment vessel the kinetic en- 
ergy must be absorbed and the containment vessel decelerated at a rate which does not 
overstress the containment vessel walls. To accomplish this, the vessel is surrounded 
by an energy-absorption material. 

Energy absorbers have been used most recently to provide soft landings for space 
vehicles on lunar or planetary surfaces (refs. 1 to 3). In this application the energy ab- 
sorbers considered were material deformation, gas bags and gas-filled collapsible 
shells. 

The requirements for a candidate energy-absorption system studied in this report 
are (1) impact speeds to 600 feet per second (183 m/sec), (2) minimum weight, and (3) 
availability of test data. An energy absorber with the potential to meet these require- 
ments is the frangible tube (ref. 4). This device was selected because (1) it has been 
tested to specific energies (amount of energy that can be absorbed by 1 lb of energy ab- 
sorber) of 38 500 foot-pounds per pound (115 500 J/ kg), which is twice that of other en- 
ergy absorbers, and (2) it has a potential to double these values. The weakest feature of 
the frangible tube is that it is essentially a unidirectional energy absorber. A single 
tube has only about ±5° angular capability. However, a low angular capability is com- 
mon to all deformable structure energy absorbers. To increase directional capability, 
these tubes must be integrated into a system capable of taking some lateral loading. 

The specific energies quoted in the literature for frangible tubes are based on tube 
weight only. In a flight system, die weight and associated support structure weight will 
add to the tube weight, thus reducing the specific energy. Only through a system analy- 
sis can these weights be defined. The purpose of this report, therefore, is to analyze 
the frangible tube, integrate it into an energy-absorption system for the protection of the 
containment vessel, calculate the total weight of the system, determine the system spe- 
cific energy, and identify which design parameters most affect system weight. The pa- 
rameters that will be studied are velocity at impact, deceleration rate, ratio of tube di- 
ameter to die radius, angular capability, number of tubes, and the specific energy of the 
tubes. 
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SYMBOLS 


2 2 

A area, in. (cm ) 

2 2 

a acceleration, ft/sec (m/sec ) 

D diameter, in. (cm) 

o 

E modulus of elasticity, psi (N/cm ) 

F force, lb (N) 

g gravitational constant, 32 ft/sec (9.85 m/sec ) 

L tube length, in. (cm) 

N number of segments 

P pressure, psi (N/cm ) 

R die forming radius, in. (cm) 

r radius, in. (cm) 

S stroke, in. (cm) 

spe specific energy, ft-lb/lb (J/kg) 

t wall thickness, in. (cm) 

v velocity, ft/sec (m/sec) 

Av velocity change, ft/sec (m/sec) 

w weight, lb (kg) 

p density, lb/in.^ (kg/m^) 

2 

ct stress, psi (N/cm ) 

Subscripts: 
b bulkhead 

d die 

ds die shank 

F f ranging 

f final velocity 

i initial inside diameter 

im impact 

o outside diameter 
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p protected package 

s shear 

sys system 

t tube 

y yield 

DESCRIPTION OF PROBLEM 

One of the phases of the crash sequence is the ground impact phase, which starts 
when the aircraft contacts the ground. If the impact occurs as the result of a midair 
collision or structural failure, the aircraft can strike the ground at any angle, thus 
making the velocity difficult to determine. Even if the aircraft is intact on contact with 
the ground, the angle of impact and the velocity can still vary over a wide range. In this 



(a) Low velocity - low angle impacts. Energy is absorbed primarily by sliding friction between plane and ground. 



CD- 10656- 10 

(b) High velocity - high angle impacts. Energy is absorbed primarily by deformation of material. 

Figure 1. - Typical types of aircraft impacts. 
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report, the containment vessel is assumed to be able to strike the ground at any angle, 
and a range of velocities is considered. 

Figure 1 illustrates the two major types of aircraft crashes: low velocity - low 
impact angle, and high velocity - high impact angle. In the low velocity - low impact 
angle impacts, the energy is absorbed primarily by sliding friction between the aircraft 
and the ground. For such an impact, the aircraft may break up into several large pieces. 
In the high velocity - high angle impacts, the aircraft kinetic energy is absorbed by de- 
forming or crushing the aircraft structure. Energy absorbers are required for this type 
of impact. 

The major characteristics of passive energy absorbers are (1) they do not have to be 
activated or triggered, (2) they are ready to absorb energy at all times, and (3) they do 
not interfere with the operation of the aircraft. The following sections explain how 
frangible-tube energy absorbers work, review experimental data, and describe the de- 
sign of an energy-absorption system. 


FRANGIBLE-TUBE ENERGY ABSORBER 
Energy-Absorption Mechanism 

Frangible tubes utilize a fragmenting tube process in which energy is absorbed 
through the force developed when a frangible tube is pressed over a die (see fig. 2). 


F 



Tube and die Initial fragmentation Final fragmentation 

CD- 10658-02 

Figure 2. - Fragmenting process. 
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(a) Tube engaging on taper of die. 



(b) Stretching of tube-strain energy due to hoop tension. 



(c) Tube curving with radius-strain energy due to bending. 


Figure 3. - Fragmenting tube process. 

CO- 10659-02 

This process has been both structurally analyzed (ref. 5) and experimentally evaluated 
(ref. 4) to determine the performance of frangible tubes. Analytically, the process can 
best be explained by the use of figure 3. There are three stages that occur in the frag- 
mentation process: 

(1) Initially the tubes engage the die on a tapered surface above the die radius as il- 
lustrated in figure 3(a). 

(2) As the tube is forced onto the die, stretching occurs and wall thicknesses are 
reduced. As a result, local failure occurs due to hoop tension in the tube (ref. fig. 3(b)). 

(3) Continued travel of the tube results in the lip of the tube rolling up within the die 
radius with subsequent fragmentation of the strips as shown in figure 3(c). 

Energy is absorbed in the fragmentation process by plastically deforming the material 
beyond its yield point. 

The amount of energy that can be absorbed is determined by the yield strength and 
ductility of the tube material and by the design of the die. Mathematically the energy 
absorbed can be expressed by 

Energy = (Force) (Stroke) = FS (1) 

Specific energy = gSl^l H Strck e ) (2) 

Weight 
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Now in the case of the fragmenting tube, a franging stress is defined as the applied 
axial force necessary for the franging process divided by the tube cross-sectional area. 
Therefore, 


°F = 


Force 


(3) 


Energy = °pA t S (4) 

and the weight of the tube f ranged is p^ A^S. Finally, the specific energy is based only 
on tube weight: 


ct f A.S ct f 

Specific energy = = — 

P t A t S p t 


Thus, the specific energy obtained from a tube is a function of the franging stress to 
density ratio of the material. With a properly designed die the franging stress ap- 
proaches the yield stress of the material. As a limit, the maximum expected possible 

specific energy is the yield stress to density ratio o /p, . For comparison of materials 

y t 

this value is often used. Table I illustrates yield stress and density values of some 

promising materials having high energy- absorption capability. Although the ct /p. 

y t 

ratios shown are encouraging in comparison with current technology, it must be noted 
that the mechanism of franging a tube is not well known. Therefore, the performance 
falls short of the potential a /p. values. 

y ^ 

Frangible tubes, however, appear to have a unique characteristic over many other 
energy absorbing devices in that the specific energies recorded in laboratory tests do 
more closely approach the maximum possible specific energy of the material a /p.. 

y ~ 

Honeycomb structures, for example, use the same materials as the frangible tube (e.g. , 
aluminum, maraging steel), but their recorded specific energies fall far short of the 
material’s o^/p potential. Measured frangible tube specific energies are, at present, 
about 60 percent of the maximum possible. 


Tube Specification and Selection 

The capability of the tube then is based primarily on the specific energy it develops. 
More specifically, it is the ratio o-p/p^. Although the franging stress is dependent on 
both the tube material and die design, with a properly designed die, the greater the yield 
strength of the tube the greater the ct f . 

7 


I 



TABLE I. - MAXIMUM ENERGY ABSORBING CAPABILITY OF CANDIDATE MATERIALS 


Material 

Yield stress, 

a y’ 2 
psi (N/cim) 

Density, 

P, 

lb/in. ^ (kg/cm^) 

Maximum energy 
absorption capability, 
°y/P, 

ft- lb/lb (j/kg) 

Maraging steel 
(18 Ni maraging) 

300 000 (206 000) 

0.289 (0.798X10 4 ) 

86 500 (259 000) 

Plastic a 

( M E M glass fabric and epoxy) 

60 000 (41 300) 

0.065 (0. 179X10 4 ) 

77 000 (230 000) 

Titanium 

(6A1-4V titanium alloy) 

126 000 (86 600) 

0. 160 (0.443X10 4 ) 

65 600 (196 400) 

Steel alloy (AISI 4130) 

210 000 (144 800) 

0.283 (0. 783X10 4 ) 

61 800 (185 000) 

Aluminum alloy (7075 -T6) 

73 000 (50 300) 

0. 100 (0.277X10 4 ) 

60 800 (182 000) 

Aluminum (2024 -T3) 

42 000 (29 000) 

0. 100 (0.277XI0 4 ) 

35 000 (104 800) 

Magnesium (AZ31B- H24) 

16 000 (11 000) 

0.0639 (0. 177X10 4 ) 

20 900 (62 500) 


a E-glass filament wound plastic material has greater yield stress. However, E-glass 
filament wound tubes were tested in ref. 9 and had lower energy absorbing capability 
than the E-glass fabric tubes. 


The candidate materials selected from table I whose franging stress would ap- 
pear to be the highest are maraging steel, titanium, and aluminum. Plastic material 
also shows good Oy values which would indicate high franging stress cr-p. 


Review of Frangible-Tube Experimental Data 

There has been some experimental work conducted by Langley Research Center on 
frangible tubes (ref. 2). Results of this effort can be summarized as follows: 

(1) The average franging stress varies directly as the ratio of the tube wall thickness 
to die forming radius t^/R. 

(2) The maximum franging stress occurs at t^/R ratios of 0.6 to 0.65 depending on 
the material. When t t /R becomes greater than these values, franging stops and the 
tube fails by splitting and buckling. 

(3) For 2024-T3 aluminum alloy tubing, the franging stress varied directly as the 
cube root of the ratio of the tube inside diameter to die forming radius D/R at a con- 
stant t t /R ratio. 

(4) AISI 4130 steel was the most efficient tubing tested on the basis of specific energy 
absorption. 
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(5) Energy-absorption capabilities of various materials increase with increasing ma- 
terial yield stress. 

For a high performance tube design, the material selection and die design play an 
important role. In reference 4, AISI 4130 steel tubing demonstrated the highest specific 
energy of the materials tested. When operation was at a t^./R ratio of 0. 6, a specific 
energy of 37 000 foot-pounds per pound (111 000 JAg) or 60 percent of its yield was re- 
corded. If the franging stress were equal to its yield, energies of about 62 000 foot- 
pounds per pound (185 000 JAg) would be possible when using AISI 4130. If maraging 
steel tubes were franged at their yield stress of 300 000 psi (206 500 N/cm ), specific 
energies to 86 500 foot-pounds per pound (259 500 JAg) would be possible. 

As cited in reference 4, however, the final consideration is the ability of the tube to 
fragment. High yield stresses do not ensure a fragmenting process. In some cases, 
tube buckling can occur prior to fragmentation. Unfortunately, the mechanical proper- 
ties of a material required for fragmentation are not well understood, and, as a conse- 
quence, the maximum Op obtained to date is about 60 percent of Oy. 


Frangible-Tube Design 

The two basic requirements in designing a frangible tube for an energy absorbing 
system are to (1) absorb the kinetic energy of the package to be protected and (2) to pro- 
vide the desired deceleration rate during this process. These are independent require- 
ments. The deceleration during slowing down is related only to the initial impact velocity 
and the allowable stopping distance; namely, 


a = 


V? 

lm 

2S 


( 6 ) 


The weight of the package does not affect the deceleration nor does the specific en- 
ergy of the tube (assuming all the kinetic energy is dissipated by the energy absorber in 
the stroke specified). The weight of tube required to decelerate the package is a function 
of the kinetic energy of the impacting package and the tube specific energy; that is, 

Tube weight (lb) = Kinetic energy (7) 

Specific energy of tube 


The kinetic energy is expressed by 


KE=±™ vjL 

2 g im 


( 8 ) 
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where w is weight of impacting package plus the weight of the tube and die. With equa- 
tions (6) and (7) satisfied, the tube and die are dimensioned to be compatible with force, 
t^/R, and cr-p requirements. In addition, a length to diameter ratio is imposed on the 
tube to ensure that it will not buckle. A detailed derivation of these relations is pre- 
sented in appendix A. 


Frangible-Tube Operating Angular Limitation 

The frangible tube is quite sensitive to the direction the force is applied to the tube 
axis (estimated angular capability, ±5°). As this angle increases, the tube will cease 
franging and will begin buckling. Therefore, if the angle of the applied force is unknown, 
several tubes must be distributed throughout the solid angle. This tube grouping (here- 
after called a segment) performs in the same manner as a single tube; that is, the total 
number of tubes absorb the total impact energy. 


ENERGY ABSORBER SYSTEM DESIGN 

Thus far the emphasis has been concentrated on the design of the frangible tube and 
die. As the amount of energy to be absorbed increases because of the high impact ve- 
locities and/or large weights, multiple tubes and multiple layers are required (due to 
L/D limitations imposed), resulting in additional superstructure for support. This 
more complex arrangement is defined as an energy-absorption system. 


System Design Prxedure 

In the case of a nuclear airplane, the protected package is a reactor containment 
vessel, spherical in shape. When designing a frangible tube for protecting this package, 
a one-tube design would result in an operating angle limitation of ±5°. It would, there- 
fore, be desirable to use multiple tubes and support the tubes laterally. With multiple 
tubes the energy absorber operating angle may conceptually be increased to ±30°. The 
requirement for protecting the containment sphere over a wide range of impact angles 
still dictates that the segments be repeated many times. An operating angle of ±30° re- 
quires 16 segments to protect the vessel in all directions. 

A requirement of low g's in the package to be protected results in long strokes (en- 
ergy absorber lengths). A limitation on the length to diameter ratio of the tube to pre- 
vent tubular column buckling is also added to the system complexity. Although no fur- 
ther tubes are necessary, a second bulkhead or structure plate is necessary for lateral 
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support at an L/D of say 10 . This adds to the system weight through additions of sup- 
port structure. 

A system design analysis, therefore, involves not only the tube and die design for 
optimum absorption of energy but the design of supporting structure for the transmission 
of the deceleration loads. All of these components add to the weight of the system and 
energy to be absorbed. Increasing the weight increases the kinetic energy and, accord- 
ing to equation (7), increases the number of tubes required. 


System Design Characteristics 


For the system analyzed in this report (i. e. , a nuclear core surrounded by shielding 
and containment vessel), the design characteristics are as follows: 


Weight of protected package, lb (kg) . . 

Deceleration, g 

Impact velocity profile, ft/sec (m/sec): 

In frontal direction 

In all other directions 

Sphere diameter, ft (m) 

Number of segments 

Tube material 

Tube specific energy, ft- lb/lb (J/kg). . 

Design ratio, t^./R 

Design ratio, D/R 

Tube L/D limitation 

Die design 

Die material 

9 

Die yield stress, psi (N/cm ) 

o 

Die shear stress, psi (N/m ) 

Bulkhead wall thickness, in. (cm) . . . 

Bulkhead material 

Impact surface 


50 000 to 400 000 (22 600 to 181 000) 
140 to 600 

250 to 600 (76 to 183) 

250 (183) 

12 (3.67) 

18 to 30 

. . Maraging steel (18 Ni maraging) 

40 000 to 90 000 (120 000 to 270 000) 

0.6 

8 to 12 

10 

See fig. 10 

Plastic with steel liner 

100 000 (68 900) 

100 000 (68 900) 

0.250 (0.635) 

Plastic 

Granite 


System Design 

In designing a system, the tube and die arrangement chosen reflect the configuration 
of minimum weight and reliable operation. 

Review of candidate system designs. - In figure 4 a multilayer 4 it protection (omni- 
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directional) design is shown. Here the tubes are mounted on the sphere first and then 
the dies. Figure 5 shows schematically other arrangements. Figure 5(a) shows the dies 
reversed to that shown in figure 4. Figure 5(b) shows tapered tubes with one die and fig- 
ure 5(c) shows an untapered tube. A detailed analysis of the directional problem may 
well show distinct advantages of one or more of these arrangements. They are pre- 
sented here as but a few alternates to the arrangement analyzed in this report. 

Design analysis . - The system design considered in this report contains the tube 
and die arrangement of figure 5(a). In this configuration the die is attached to the sphere 
with the tube extending radially outward. 

The tubes are then attached for lateral support to a bulkhead or structural plate. 
Where multiple layers are necessary, a die is mounted on the bulkhead with a second 
tube extending radially outward to a second bulkhead. 

The analysis has been programmed for a IBM 7094 computer (see appendixes A 
and B). This analysis provides for multiple grouping of tubes, each of which contains 
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multiple layers where required. Stress calculations are provided for the dies of each 
tube. Estimated weights are calculated for the additional superstructure that is re- 
quired. 


PARAMETRIC ANALYSIS OF SYSTEM 

In designing a frangible-tube system there are design parameters that can be ad- 
justed for the optimization of a minimum weight system. The purpose of this parametric 
analysis, then, is to determine (1) the minimum weight possible for the system, (2) the 
design parameters having the greatest effect on this weight, and (3) the design parame- 
ters having a lesser effect on this weight. 


Description of Parameters 

The parameters that form a part of the system to be analyzed are the number of 
segments, number of tubes per segment, specific energy of the tube, impact velocity, 
deceleration during impact, and ratio of frangible-tube diameter to die radius. 

Number of segments . - Each group of tubes mounted radially on the sphere capable 
of absorbing the energy of impact is referred to as a segment. The greater the number 
of segments the greater the system weight but the smaller the included angle within 
which the tubes are required to operate (see fig. 6). 

Number of tubes . - Within a segment there can be a variation in the number of 
tubes. With a few tubes their diameters, wall thicknesses, and dies are large. With 
many tubes their diameters, wall thicknesses, and dies are small. 

Specific energy of tube . - The specific energy of the tubes dictates the total weight 
of tubes of each segment because the product of the total weight of the franged material 
and its specific energy must equal the kinetic energy the segment is to absorb. (The 
specific energy of the tube considers only the tube weight - not the die weight. ) 

Impact velocity . - For a given protected package weight the kinetic energy to be ab- 
sorbed increases as the square of the impact velocity, thus, becoming an increasingly 
major contributor to the energy absorber weight. 

Deceleration during impact . - At impact the impact velocity is reduced to zero at a 
predetermined rate or deceleration g. The higher the g values, the smaller the stroke. 
High g values, however, result in high forces, heavy support structures, and thick con- 
tainment vessels. Low g's result in long strokes, resulting in more layers of energy 
absorbers (due to L/D limitations). 

Ratio of frangible tube inside diameter to die forming radius, D/R. - The variation 
of this parameter affects primarily the die weight. As the tube diameter increases the 


14 




Figure 6. - Energy absorber weight and segment in- 
cluded angle as function of number of segments. 
Payload, 200 000 pounds (90 600 kg); deceleration, 
300 g' s; velocity profile, 400 feet per second 
(122 m/sec) +30° from front and 250 feet per sec- 
ond (76 m/sec) in all other directions. 


die diameter also increases, resulting in a weight addition as a function of the die diam- 
eter squared. Upon examination of all of these parameters a minimum weight system 
can be obtained. 


Design Constraints and Assumptions 

A minimum-weight energy-absorption system has been designed for protection of a 
typical nuclear reactor power source. The reactor would generate 300 megawatts of 
power and be applicable to a 1 million pound (4.5x10^ kg) aircraft. The reactor core and 
shielding are surrounded by a large spherical containment vessel 12 feet (3. 67 m) in di- 
ameter. The following assumptions were made in the design analysis: 

(1) The frangible tubes are mounted on the surface of the sphere extending radially 
outward with the tube and die arrangement of figure 5(a). 

(2) The energy- absorption system possesses an omnidirectional capability. This 
condition requires that the ’ 'protected package” be surrounded by energy absorbers in a 
47r containment. 


15 



(3) The maximum impact velocities occur only in the frontal direction (±30° solid 
angle). The remaining directions are capable of absorbing energies at impact velocities 
no greater than 250 feet per second (76 m/sec). The aircraft cruises at a velocity of 700 
to 800 feet per second (214 to 244 m/sec). 

(4) The segment impacting absorbs only 75 percent of the total energy of impact. 

The remaining 25 percent is absorbed by some of the redundant segments. During im- 
pact, the franging of the tubes in the impacted segment also results in the folding or col- 
lapsing of adjacent segments. The 75 to 25 percent assumption takes this energy absorb- 
ing process into consideration. 

(5) No energy is absorbed by the air frame. 

(6) No energy is absorbed by the containment vessel. 

(7) Frangible-tube dies are fabricated from plastic material with a metal- lined 
forming radius. 

(8) Plastic bulkheads are provided to absorb the lateral loads. Where strokes ex- 
ceed an L/D of 10, multiple bulkheads are provided. 


RESULTS 

The preliminary analysis that was conducted is divided into two categories for dis- 
cussion: (1) those design parameters having the greatest effect on the system weight, 
and (2) those design parameters having the lesser effect on system weight. All parame- 
ters studied, regardless of the category they fall in, were varied over a range consid- 
ered to be within a reasonable extrapolation of existing state-of-the-art technology. The 
ranges of these variables were included in the System Design Specifications and Require- 
ments section. 


Design Parameters Having Greatest Effect on System Weight 

The system design parameters which have the greatest effect on the energy absorber 
system weight are the number of segments (and subsequently the angular dependence of 
each segment), the specific energy of the frangible tubes, the frontal impact velocity, 
and the deceleration of the protected package. 

Number of segments . - Figure 6 contains a plot of the energy absorber system 
weight in percent of protected package weight against the number of segments in the sys- 
tem. Specific energies of 90 000 foot-pounds per pound (270 000 J/kg) (a theoretical 
maximum value using maraging steel tubes) and 60 000 foot-pounds per second (180 000 
J/kg) were considered. In addition, the segment included angle is plotted as a function 
of the number of segments. This angle represents the solid angular protection that one 
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segment provides for the sphere. It also represents the maximum angular load that is 
applied to tubes in that segment. 

The number of segments were arbitrarily varied over a large range. The result 
was that the energy absorber system weight increases sharply with the increase in the 
number of segments especially at lower specific energies. The segment included angle, 
however, decreases with increasing number of segments. The need for reducing lateral 
loading on the tube within the segment by increasing the number of segments thus is off- 
set significantly by an increase in system weight. 



60 180 300xl0 3 

Specific energy, J /kg 

Figure 7. - Energy absorber weight as function 
of tube specific energy. Payload, ZOO 000 
pounds (90 600 kg) ; deceleration, 150 g T s ; 
velocity profile, 250 to 400 feet per second 
(76 to 122 m/sec); segment operating angle, 
46°; number of segments, 24; tube diameter 
to die forming radius ratio, 10. 


Frangible-tube specific energy . - The effects of varying the frangible-tube specific 
energy in a system of 24 segments and at a 400 foot per second (122 m/sec) frontal im- 
pact velocity protection are illustrated in figure 7. At a frangible- tube specific energy 
of 90 000 foot-pounds per pound (270 000 JAg) the energy absorber system weight is 
46 percent of the protected package weight. This weight increases to 76 percent with a 
one-third reduction in specific energy to 60 000 foot-pounds per pound (180 000 JAg). 
Further reductions in tube specific energies result in prohibitive weight penalties . The 
need for high specific energies is well illustrated in figure 7. 

Frontal impact velocity . - The effects of increasing the frontal velocity under the 
same conditions of frangible-tube specific energy and 24 segments is illustrated in fig- 


I 
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90 120 150 180 

Frontal impact velocity, m/sec 

Figure 8. - Energy absorber weight as function of frontal impact ve- 
locity. All remaining directions are designed for impact velxity 
of 250 feet per second (76 m/sec). Payload, 200 000 pounds 
(90 600 kg); deceleration, 300 g's ; segment operating angle, 46°; 
number of segments, 24; number of tubes, 40; specific energy, 

90 000 foot-pounds per pound (270 000 J/kg). 


ure 8. Increasing the frontal impact velocity from the minimum value of 250 to 600 feet 
per second (76 to 183 m/sec) increases the system weight by 68 percent. 

Deceleration of protected package . - For the deceleration parameter, a minimum 
system weight is obtained at 250 g’s. This can best be explained as follows. With a re- 
quirement of low g’s the stroke (length of tube to be franged) must be large, thus re- 
quiring multiple layers. These layers add more dies and bulkheads to the system and 
thus more weight. Now as the deceleration values increase, the stroke is reduced. 

This condition eliminates bulkheads and dies, thus reducing the system weight. The in- 
crease in deceleration values, however, increases the force per tube requiring larger 
tube diameters, wall thicknesses, and die diameters, which thus increases the die 
weight. For example, a typical die weight at 140 g's of deceleration is 5.5 pounds 
(2.5 kg); this weight increases to 50 pounds (22.6 kg) at 600 g's. This increased die 
weight offsets the weight reduction due to fewer bulkheads. 


Design Parameters Having Lesser Effect on System Weight 

The second area of investigation was the design parameters which had the lesser 
effect on the energy absorber system weight. Two variables were analyzed - namely, 
the number of tubes per segment and the frangible-tube diameter to die radius ratio. 
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Figure 9. - Energy absorber weight as func- 
tion of number of tubes per segment and 
package deceleration. Payload, 200 000 
pounds (90 600 kg); velocity profile, 250 to 
400 feet per second (76 to 122 m/sec) ; seg- 
ment operating angle, 56°; number of seg- 
ments, 18; specific energy, 90 000 foot- 
pounds per pound (270 000 J/kg). 


Number of tubes per segment . - In the case of the number of tubes per segment, 
illustrated in figure 9, the curve is fairly flat showing that the minor weight penalty for 
increasing or decreasing the number of tubes per segment would be dictated by geometry 
and mechanical design considerations. 

Frangible-tube diameter to die radius ratio . - The ratio D/R was varied between 
8 and 12 at a constant t/R ratio of 0.6. These D/R values were recommended in ref- 
erence 4 for best franging performance. The calculations showed that this parameter 
had no effect on the system weight. It was expected that as this ratio increased a reduc- 
tion in system weight would occur, since large D/R ratios result in larger tube diame- 
ter, longer length tubes for a L/D of 10, and, subsequently, fewer layers, dies, etc. 
Between the range selected (8 to 12) this effect was either too slight to notice or was off- 
set by increasing die weights. 


DISCUSSION OF RESULTS 

The parametric analysis reveals that using conventional energy absorbers, such as 
frangible tubes, results in high system weights . As an example, a typical 300-megawatt 
powerplant system design being studied at Lewis for application to a 1 million pound 
(4.53x10^ kg) nuclear airplane has the following specifications: 
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Weight of powerplant (containment vessel, shielding, and core), 


lb (kg) 200 000 (91 000) 

Impact deceleration, g 300 

Impact velocity profile, ft/sec (m/sec) 

In frontal direction 400 (122) 

In all other directions 250 (76) 

Containment vessel diameter, ft (m) 12 (3.66) 

Number of segments 24 

Segment operating angle, deg 46 

Number of tubes per segment 40 

Specific energy of tubes, ft- lb/lb (J/kg) 90 000 (270 000) 


The dies, bulkheads, etc. of the system were assumed to be fabricated from a high 
strength plastic material. 

From the parametric curves, the resultant weight of this system is approximately 
50 percent of the powerplant. The specific energy of the overall system is 1800 foot- 
pounds per second or only 2 percent of the assumed tube specific energy. The principal 
contributors to these large weights are the segment redundancy and the required super- 
structure (dies, bulkheads) which offset the high specific energy obtained from the tubes. 

A more detailed design analysis of a frangible-tube system was not conducted. The 
reason is that, although the parametric analysis conducted herein was preliminary in 
nature, the design assumptions made were optimistic rather than conservative. As an 
example, tube specific energies to 90 000 foot-pounds per pound (270 000 J/kg) were 
used, plastic dies were assumed which are currently undeveloped, and thin plastic bulk- 
heads were utilized for the lateral tube support and multiple layer application. A more 
detailed structural analysis may tend to increase the weight of these components and 
subsequently the overall weight of the system. 

One area of consideration for weight reduction is the utilization of other reactor and 
airplane systems components as energy absorbers. An example would be the integration 
of the shielding as both a shield and energy absorber and use of the airframe to absorb 
some energy. 


CONCLUDING REMARKS 

The nuclear airplane must be safe. To meet this safety requirement the nuclear 
power source is protected at impact against release of fission products. Frangible-tube 
energy absorbers are candidate devices for this protection. Their recorded high specific 
energies of 38 500 foot-pounds per pound (115 500 J/kg) together with the potential for 
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even higher specific energies are their major advantage. This specific energy, how- 
ever, reflects only the tube weight and does not include die weight and associated struc- 
tural weights. 

Therefore, a parametric study was conducted for determining the weight of one 
energy-absorption system concept utilizing frangible tubes. The system typified that 
which might be applied to a nuclear airplane. The results of this study are as follows: 

1. Decreasing the number of segments surrounding the sphere or increasing the 
specific energy of a frangible tube substantially decreases the energy absorber system 
weight. Variables of lesser effect but still of significance are the frontal impact velocity 
and the deceleration of the protected package. 

2. Changing the number of tubes per segment or changing the design ratio of tube 
diameter to die radius has a minor effect on system weight. 

The parametric analysis revealed that with frangible tubes the system weight in all 
cases was large. A typical 300-megawatt powerplant for a 1 million pound (4. 53x10^ kg) 
nuclear airplane would require an energy-absorption system weighing 50 percent of the 
powerplant weight. 

Based on the parametric analysis, the following recommendations are made: 

1. Other candidate systems should be studied before frangible -tube systems are 
looked at in more detail. The reason is that although the parametric analysis conducted 
was preliminary, the design assumptions made were optimistic rather than conservative. 
That is, the performance data of the tubes and support structure were assumed at their 
theoretical maximum values, resulting in a minimum system weight. 

2. Methods must be studied in which the energy absorber is integrated with the 
shield, structural material, etc. In this manner the system components serve both as 
they are intended to be used during normal system operation and as an energy absorber 
during impact. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 17, 1969, 

126-15. 
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APPENDIX A 


DESIGN ANALYSIS 

The design analysis developed for the energy-absorption system calculates (1) the 
size of the tube necessary to absorb the kinetic energy of the protected package and (2) 
the weight of the tube, die, and bulkhead which comprise the energy-absorption system. 
Since the tube, die, and bulkhead weights also become part of the weight of the package 
to be stopped, both calculations are dependent on each other. 


Calculating Tube Size 

To arrest a package in motion at a controlled magnitude of deceleration, it is neces- 
sary to oppose the motion with a force. This relation is 

w„ 

F = — H a (Al) 

g 

At the initial time at which the force is applied, the total energy the package possesses 
is expressed by 


1 w g 

Kinetic energy = E Vf 

2 g 1 


(A2) 


In the simple case where one franging tube provides the resisting force, the length of the 
tube franged must be 


S 


2 g F 


(A3) 


when a constant force is assumed over the distance S. 

With the equations of energy satisfied, it is now necessary to design a frangible tube 
which will provide a force F through a distance S. The franging stress of a tube is 
defined as 


a F 



(A4) 
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Some empirical relations of franging stress as a function of tube and die design are as 
follows: 


ct f = 1000(284.0 — - 58.01 

V R 


(A5) 


°F = 


D t 

1900 -A - 0.333 
R 

t t 

0.7 - -i 
R 


(A6) 


o-p = 1000 ( 2.21 — - 1. 1 1 

R 


(A7) 


Equation (A6) gives the franging stress for the 2024- T3 aluminum alloy which is taken 
from reference 4. Equations (A5) and (A7) are franging stress approximations for 
AISI 4130 steel and A231B magnesium alloy, respectively. Both equations were derived 
from graphs in reference 4. 

When the franging stress (obtained from typical eqs. (A5) to (A7) or assumed) is 
known, the tube area is obtained by rearranging equation (A4): 

A t = — (A8) 


The tube area is also related to the inside and outside tube diameters by 



Substituting (D i t + 2t^) for D q ^ and solving for the wall thickness give 


(A9) 



(A10) 


where D i ^/t^. is obtained from input parameters t^/R and Dj t /R. 

In cases where the stroke S exceeds a tube length to diameter ratio of 10, a second 
tube is used. This requires the calculation of a final velocity for the first tube. Using 
the time motion equation the final velocity is 
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where 


V 2 - = 2aS = A(V 2 ) 


(All) 


A(V 2 ) = 2aS 


(A 12) 



(A13) 


Although Vj represents the final velocity for the first tube, it is the initial velocity for 
the second tube. This calculation continues until Vj is 0. At this point the total kinetic 
energy of the package has been expended by a franging force F through a stroke S. 

An important parameter used in evaluating a frangible tube is the specific energy. 
The specific energy is defined as the amount of energy the tube is capable of absorbing 
per pound of tube. It is derived simply as follows: 


Specific energy = Energy absorbed 
Weight of tube 

For a tube force F acting through a franging stroke S 


spe = 


FS 

P t SA t 


Substituting equation (A8) for F 



(A14) 


(A15) 


(A16) 


Calculating System Weight 

The system weight consists of the tubes, dies, and bulkheads of all segments sur- 
rounding the containment vessel. The weight is calculated by sizing the tubes and dies 
in each segment as if that segment were the one in which the impact occurred. 

The weight of the tube is calculated by using the relation of the specific energy and 
the kinetic energy: 
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Tube weight = -£®- (1 + (A17) 

spe \ S/ 


Tube weight = (l + — ] (A18) 

2g spe \ S / 

The ratio 6R/S represents the percent of the tube unfranged due to the stroke bot- 
toming on the die shank (see fig. 10). The quantity 6R/S adds to the franged tube 
weight which will remain unfranged. 



Figure 10. - Frangible tube and die design dimensions. 


When calculating the weight of the die, it is necessary to determine the wall thick- 
ness of the shank t , and the shear thickness t supporting the die forming radius R 
(see fig. 10). To determine these thicknesses accurately requires a complex stress 
analysis of the tube which is beyond the scope of this investigation. Therefore, a simple 
relation to porportion the surface pressure P applied to the die shank for all values of 
tube radius, wall thickness, and yield stresses was assumed to be 


P = 



(A19) 


where P is the boundary pressure on the die shank attempting to collapse it and 


°y 18 
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the tube yield stress. The allowable boundary pressure for collapsing the shank is ex- 
pressed by 


P 




(A20) 


(ref. 8) where o-y is the die yield stress. Setting equation (A19) equal to equation (A20) 
and solving for t^g yield the cubic equation 


X 


3 



where 


(A21) 


X = 


ds 


ds 

4cr' 


C l = 


E 


ds 


The required die shank wall thickness is given by r^ = r^ g + t^ g , which, after rewriting, 
becomes 


ds 


(X+ 1) 


The shear thickness t is obtained in a more direct manner: 


*\ 


CT s 


t = 


2j r i _ l u s 


(A22) 


(A23) 
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Finally, assuming a t^/R of 0. 6 the die forming radius R can be determined 
based on the tube thickness. The remaining dimensions used for the die are shown in 
figure 10 . When these dimensions are used, the equation for computing the die weight 
is (see fig. 10) 


Die weight = Weight of die shank + Weight of die base 



( D i 2 ,t- D Ms) (,R+ V + jPc 



+ 4R + 4t t ) 2 



(R + t g ) 


- 7T 2 R 2 (D, * + 2R) (A24) 

9 1 


The last component of the system for weight estimates is the bulkhead. This equa- 
tion is simply 


Bulkhead weight = 


4?rr bVb 


N 


(A25) 


The total segment weight becomes 

Segment weight = Tube weight + Die weight + Bulkhead weight (A26) 

The total system weight is 

N 

System weight = ^ Segment weight + Protected package weight (A27) 

i=l 
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APPENDIX B 


PROGRAM LISTING 

In designing an energy- absorption system for protecting a package, the process is 
complicated in that as energy absorbers are used to absorb kinetic energy, they too add 
to the overall weight of the system. For maximum design flexibility, the analysis is 
best accomplished on the computer using an iterative process. 

The method of calculating an energy absorber using the equations of appendix A is 
summarized as follows : 

(1) Calculate franging stress (see eqs. (A5) to (A7)) or designate a franging stress. 

(2) Calculate specific energy (see eq. (A16)). 

(3) Assume a system weight. 

(4) Begin calculating segment weights. 

(5) Begin calculating the first section weight of the first segment. 

(6) Calculate the tube force (see eq. (Al)). 

(7) Calculate the tube area, wall thickness, and outside diameter (see eqs. (A8) and 
(A10)). 

(8) Calculate franging stroke when the tube L/D = 10. 

(9) Calculate the final velocity of the layer (see eq. (A13)). 

(10) Calculate the tube weights, due weights, and bulkhead weight of the section (see 
eqs. (A18), (A24), and (A25)). 

(11) If the final velocity calculated in step (9) is not zero, return to step (5) and cal- 
culate the next section. Repeat until final velocity is zero. 

(12) Sum all section weights for the segment (see eq. (A26)). 

(13) Return to step (4) and calculate next segment. Repeat until all segments have 
been calculated. 

(14) Sum weights of all segments for system weight (see eq. (A27)). 

(15) Use new system weight for step (3) and repeat steps (4) to (14) until convergence. 
Convergence occurs when two successive system weights are within convergence toler- 
ance. 

The computer listing of the above calculations follows. All input and output are in 
U.S. customary units. 
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t P AR AMF TR I C CCDF SYSTEM 2 

C 

n IWFNS ICN DOT! 3C, 10) ,D I ( 30,10) ,VELF (30,10) ,F (30,10) ,T( 30,10 ) ,XL( 30 
I. 10).G( 30, 10 ), XNOT( 30 , 10 ) ,R ( 30 , 10) , S ( 30 , 10) , TW ( 30 , 10 ) * SWT ( 30 , 10 ) , 
20 IFWT( 30, 10 ) ,BIJWT( 30, 10) ,T$( 30 ,10) ,TP(30 ,10) , DIF (30, 10 ) 

DIMENSION VELX( 30 ) , SYSWT ( 30 ) , PRE V ( 30) ,N0(3C) 

COMMON K, J, STGF.SIGO.T.DI , S IGS ,APT ,ED , PI ,TP , CH ECK P , T S 
C 

READ (5,101) K S IGMA * NEIL ,N0 S 

REAR (5, 401) ( IXNDTI J,K) ,K=1,N0L) . J = 1 ,NOS) 

PEAR (5,001) ((G<J,K),K-1.N0L) ,J=1 ,NOS) 

READ (5,901) ( VELX ( J ),J=1,NCS) 

R F AD (5,901) R HO l , RH02 ,RH0'3 , BH T 
10 1 FflRMAT (BE 10) 

901 FORMAT (3F10.7) 

907 FORMAT (5020.9) 

904 FORMAT (11H1SFCTI0N Nt),10X,8H LAYER N0,12X,11HN0 OF TUBES ,9X , 12HACC 
1FLFRAT ION ) 

WRITE ( 6,900 
DO 903 J= 1 , NCS 
no 9 C9 K= 1 , NEL 

909 WRITF ( 6, 907 ) J , K , XNO T ( J , K ) ,G ( J , K ) 

90 B CCNTTNUF 

980 FORMAT (HHJSECTIUN NO , 10X , 1 1H INI T I AL VEL) 

WRITF (6, SBC 1 

nr SF.i j= i , n c s 

9 B 1 WRITF (6, 907) J.YELX(J) 

510 READ (5.901) 00 VR , TO VR , PL . S IG , SI G 0 , S I G S , ED 

910 FORMAT (BH.JRFOl = G 1 5 . 5 , 3 X , 7HRH 02 = G 1 5. 5 , 3 X , 7 HRH03 = GI5.5,3X,7HD 
lOVR = G15.5/EHJTOVR = G 1 5 . 5 , 3X ,1 OH PA Y LOAD =G 1 5 . 5/ 1 2H JFRAN STR = G 
21 5.5.3X. 12HDIE YIELD = G 1 5 . 5 , 3 X , 12 HDI E SHEAR = G 1 5 . 5/ 1 5 HJD I E MODUL 
3UR = G15.5) 

WRITF (6,910) RH01, RHO ? , RHO 3 , D 0 VR , TDVR , PL , SIG , SIGD ,S IGS »ED 
C 

,10 = 0 
PRF = 0.0 
P 1= 3. 161 69 

DO 7 1 = 1, NOS 
7 SYSWTI I ) = PL 
7 10 CCNT lAItJF 
NX! = 0 
S W T T = 0.0 

or 7ii j = i , no s 
nr 7ii k = l.NOi. 

711 S { J « K ) = 0. 

c 

r CAICUI.ATF FR AN G TNG STRFSS 

C 

GO TO ( 71, 72,73,74,95) , K S IGMA 
71 S 1GF= ( 284 .0*TOV/R - 58 . C 1 * 1 . OE + 03 
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C,C TO 9 5 

7 2 S TGF = 1 .9*1 .CF + C3*00VR**.33?/< .7-TCVR) 
Gf TO 95 

7 3 STGF = ( 2.21HHVR-1, 1 )*1.0E + 03 
GO TO 95 

74 SIGF = SIG 
or TP 95 

9 5 CrNTINlJF 
C 

C CALCULATE SPECIFIC ENFRGY 

C 

SPF= SIGF/(RFni*12„C) 

0 F PUG S IG F ♦ S PF 
MX= 0 

0 *** *** # « * # # 

C SFCTTON ITFRATICN LOOP 

C 

C 

OC 7PO J= 1 ,NCS 
C 

N X= ] +N X 
D F BUG NX 
DFPUf- SYSU T (J ) 

C 

NCU ) = 0 
NK= 0 

C sCs****^***#* **** * * ^ £ £ * *#2*£ 

f LAYFR ITERATION LOOP 

C ****** t ******* x-**#*#* **#**% #%%****** 

nr 5 CO K = 1 • N GL 

r 

IF (J-l) 4,4,5 

4 IF ( K- II 2,2,3 

3 SYSWT(J) = SYSWT(J) - SUT(J f K-l ) 

I F( SYSWT( J ) .LT.PL) SYShT(J)= PL 

2 GC TP 8 

5 Cf NT INUF 

IF (K-l ) 3 , F , J 0 

10 SYSWTfJ) = SYSWT(J) - ShTIJ, K-l ) 

C 

C CAI rui A TF FORCE 

3 F ( J , K ) = S Y Sift 7 ( J ) * G(J*K) 

C 

f CALCULATE TUFF AREA 

A P "F - F ( J , K ) / (SIGF * XNCT(JfK)) 

OFPUG F(JtK) , A P T 
C 

f CALCULATE TUBE THICKNESS 

DC VT = nnVR/TMVR 

T { J » K ) = SQ*T(AP1/(PI*(D0VT + 1.0))) 

r 

C CALCULATE INSTOF DIAMFTER 

r 

n i ( j ,k ) = cnv i* t ( j , k ) 

R (JtK)= T ( J t K ) / TO VR 
D F BUG OrVT,T (J,K)»DI (J.K) 
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^ ^ n n r> r> n r> n r> n 


calcui atf outside diameter 
' nrr(j,K)= oilj.k) + 2 .c*t(j,kj 

M C ( J ) = 1 + NO ( J ) 

12 XLOVD= 10.0 

XL( J,K )= XLDVD * D0T( J,K) 

S (J,K)= XL ( J,K > /12.0 
r> ELV = 2. * G(J,K) * SIJ.K) * 32.2 
DEBUG S ( J , K ) ,0P T(J»K)»DELV 

CALCULATE p INAL VELOCITY 

I F ( K — 1 ) 15, 15, 14 

15 V FL FX= VPLX(J)**2 - OELV 
I F( VELFX.LT .C.O) SO TO 17 
VFLF(J,K)= SCRT(VELFX) 

D FPtJG VPI F( J ,K ) 

SC T n 400 

17 D FL V= VEL X ( J ) **2 
NK= 1 

S(J,K)= OELV /( 2 . * G(J,K) * 32.2) 

V FL F ( ,1 , K ) = C .0 
SC TP 400 

14 V FI FX = VFI E ( J , X-l 1**2 - OELV 
I F ( V F| PX.LT.f.O) GO TO 18 
VFI c(.|,k)= SCRTt VELFX) 

DEBUG VEL F( J ,K 1 
GC TP 400 

IB 0 FL V= V FL F ( J , K- 1 )**2 
M K = 1 

S ( J . K > = OELV/f B.*G(J,K)*32.2) 

V EL F{ J ,K ) = C .0 
DEBUG SIJ ,K ) 

C Al GUI 4TB TU P c k T 

40 0 TF(K-l) 2 C , 2 C , 2 1 

20 TW(J.K) = t By SUT( .) )* (VFLX( J>**2 - VFLF ( J * K ) ** 2 ) ) / ( 64 . 4 *S PE ) * ( 1.+6 
1. *P ( J. K )/( S< J. K 1*12. ) ) 

DEBUG TW( I ,K 1 
GC TD 20 

21 TW(J,K) = ( SYSWTI J ) * ( VEl F ( J,K- 1 >**2-VELF ( J ,K)**2) ) / (64 ,4*SPE)* 

)( l.+6.*R(J»K )/I S(J,K)*12.) ) 

DEBUG T W ( J , K ) 

20 GCNT IMIIE 

C. AIGUL ATF DIF WT 

CALI DIPSTR 
C 

2 5 0 IFW=RHC2*. ?E*P I*< (01 < J,K )**2-DID< J,K )**2) *(7.*P ( J ,K)+TS ( J . K) )+{ {< 
ID T< .J,K. )+4.*R U , K)+4.*T( J,K))**2-DI < J,K)**2)*<R( J ,K)+TS { J,K) H-2.0* 
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r> r> n nrjr>r>nr>nr> n n 


?P I*R{ J.K )**?*( n I ( J ♦ K )+2.*R( J.K) ) ) 
DFRUG DIEW 

DIEWT(.I,K) = CIEW * XNOT ( J * K ) 

CALCULATF BULKHEAD WT 
RAOX = 0. 

nr aob i =i.kol 

405 RAOX = RADX + S( J. I ) 

ARF AX = 4 . / FLOAT { NO S ) *PI*(RAQX + 
RUWT(J.K) = AREAX*BHT*PH03 *144. 
DEBUG BUWT(J.K) 

CALCULATE SYSTEM WT 

SWTIJ.K )= Tto(J.K) + DIEWT(J.K) + 
DFPIJG SWT(J.K) 

31 IF (NK) 5CC, 500.600 

500 CONTINUE 


600 CENT INU F 

T UW= 0.0 

niFT= 0.0 

RUT = 0.0 

NCX = NOt J 1 

on 30 K = 1 . NO X 

T UW = TUW + TW(J.K) 

0 ] FT® OIFT + DIFWTt J.K) 

30 RUT® BUT + RLWT(J.K) 

I F ( J — 1 ) 40.4C.41 

40 S YSWTI J )® ( 0 IF T + BUT 4 - TUW) + PL 
GO TO 42 

41 SYSWTl J )= (DIET + BUT + TUW) 

4? DFRUG SYSVii T ( J) 


700 CONTTNUF 

DO 701 1 = l.NOS 

701 SWFT = SWTT •* SYSWT ( I ) 

frr= ARS(SWT) - °R E ) /SWTT 
PPF = SWTT 

IF (FRR.GT..C1) GO TO 702 
GO TO 705 
70? NXL® 1 
70 3 CONTINUE 

0 F PUG SWTT 

nr 7fii i = i.nos 

781 S YSWTt I 1 = SWTT 

I F( SWTT .GT.5.0F+C5) GO TO 7C5 

jr - ! + jo 

IF (NXL) 7C5.705.71C 


6 . 01**2 


BUW/TI J.K) 
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705 CONTINUE 

912 FORMAT ( UHLSFCTION NO . 1 OX ♦ HHL A YER NO , 12 X , 1 4HTUBF T H I C KNESS , 6X f 1 2HD 
lUTSlOF HI AM » 3X * 9 HF INAL VEL) 

WPITFf 6, 9 1 2 ) 

914 FORMAT ( 6G2C .9 ) 

nr 510 j= i * ncs 

NOX = N0( J ) 

nn 9io k= l • ncx 

910 WRITE (6,914) J , K , T ( J , K ) . DO T ( J . K ) f VE L F ( J , K ) 

921 FORMAT! 1 1 HI. S FC T ION NO . 10 X , 8HLA YER NO , 12 X , 10H0 1 E R AD I US . 10 X , 5 HF CPCE 
1.15X.11HTUBF L ENGTH , 1 8 X , 1H S ) 

WRITE ( 6, 9? 1 ) 

00 922 J - l * N f S 
NCX = NO ( J ) 
nn 922 K= 1 • NCX 

922 WRITE (6*914) J • K f R ( J * K ) , F ( J , K ) , XL ( J , K ) , S < J , K ) 

801 FORMAT (UHLSFCTION NO , 10 X * 8HLA YE R NO . L 2 X , 10H S E CT I C N WT,lOX t 6HDIE 
1WT, 14X, 1 1HBULKHEAD WT,9X,7HTU8F WT) 

WR ITF ( 6. 8C 1) 

OP 802 J=1,NDS 
NOX = N0( J ) 
nc 807 K= 1 f Nl'JX 

802 WRITF (6.914) J , K , S WT ( J , K ) , D I E k T ( J • K ) , BU U T ( J , K ) , TW ( J , K ) 

815 FORMAT (UHLSFCTION NO , 1 0 X . 8HL A YF P NO . 1 2 X * 14H 0 I F WALL THICK.6X.L5H 
1 D IF SHF AR TH IC K ) 

WRTTF ( 6. 815 ) 

Of 816 J= 1 . NOS 
NCX = NO ( J ) 
nr 816 K=l. NCX 

816 WRITF (6.914) J . K . T P ( J . K ) . T S ( J . K > 

800 FORMAT (17HJT0TAL SYS WT = G 1 5. 5 , 4 X . 1 6HF R ANGE STRESS = G15*5,4X,1 
18F SPECIFIC ENERGY = G15.5) 

WRITF (6.808) S W T T . S I G F , SPF 
WRTTF ( 6, 9] 4 ) JO 
G C TO 5 10 


SIBROUTINF CIESTR 
C 

c 

niMFNSION DOTf 30,10 ),DH 30,10) ,V£LF( 30,10 ) ,F( 3D,10 ),T( 30,10 ),XL( 30 
)• 10 ),G( 30, 1C ), XNOT ( 30 f 10) ,R<30 ,10) «S (30 , 10) ,TW (30 , 10) ,SWT( 30,10), 
POTFWK 30, 10) ,BtJUT( 30,10 ,TS( 30.10) ,TP( 30,10) ,010(30,10) 

RTMFNSIPN VFLXt 30) , SYSViT(3C),PREV(30) ,N0(30) 

COMMON K, J,S IGF, S I GO, T, 01 ,SIGS , APT, ED, PI , T P , CH ECK P , T S 
C 

X C 1 = A.*STGO/FD 
DEBUG XC1 
A = 1.0/XC1 
DEBUG A 

Xf? = -SIGC/< 4**T(J ,K)/ni ( J,K)*SIGF*XC1 ) 

DEBUG XC2 
B= XC? 

nUM= SORT ( B**2/4.+A**3/27. ) 

DERUG OUM 

B IGA = (-B/2.*DUM)**.233 
OF RIJG B IGA 
V A R = (-B/2.-CUM) 

B IGR = SIGN(ARS( VAR )**.333,VAR) 

DEBUG BIGB 
X = BIGA + BIGB 

C E F C K P = 1 ./X*SICD/( l. + XCl*X**2> 

DEBUG CF^CKP 
C 

TP(J,K)= OT ( J»K ) /( 2.*( X+ 1. ) ) 

DEBUG TP ( J f K ) 

nrn( j,k )= oi (J * k>-?.*tp( j.k) 

c 

C CALCULATE DTE SHEAR THICKNESS TS 

C 

X l H A0= SIGF*APT 

TS( J,K)= XLOAD/I SIGS*P 1*01 ( J,K) ) 

DEPUG TS( J,K ) 

C 

R FTIJPN 
END 
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